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1. Introduction 

The lattice simulation with the overlap fermions ^ keeps chiral symmetry for any number 
of dynamical flavors. It implies that the chiral properties of the numerical data are described by 
the continuum chiral perturbation theory (ChPT). It therefore makes sense to test the consistency 
between ChPT and QCD. With this motivation, we performed numerical simulation using overlap 
fermion action with two dynamical flavors [g]. In particular, we carried out the calculation of light 
meson spectrum. We fitted the data to the next-to-leading order (NLO) formulae with different ex- 
pansion parameters m q , nr % and = m 2 n / (Anfx) 2 , all of which should give equivalent prescription 
at NLO. As a result, we found that the fit curves start to deviate around the scale of kaon mass 
niK- Another important observation was that only the § -fit reasonably describes the data beyond 
the fitted region. We also demonstrated that next-to-next-to leading (NNLO) effect is necessary to 
describe the data around and beyond m^. 

Based on these findings, in this article, we present an extension of this study to Nf = 2 + 1. 
Using overlap fermion action with 2+1 dynamical flavors, we calculate m 2 n /m U( i, m\/m s d, f% and 
fx, where m u d is the u-d degenerate mass and m s d = \{m s + m^). Since these results depend on 
strange quark mass, the chiral extrapolation should be performed with a fit ansatz valid beyond the 
scale of kaon mass. The only possibility is the fit with the NNLO ChPT effect taken into account. 

After explaining how we obtain the data points briefly in Section we present in Section || 
the chiral extrapolation with NNLO ChPT formulae to treat both pion and kaon sectors on an equal 
footing. Preliminary results of physical values are also given in this section. As a consistency check 
of our analysis, in Section |], we perform the chiral extrapolation with a reduced theory in which 
the strange quark is integrated out. Section ||] contains a brief discussion about what to be done to 
get to final results on this work. 

2. Data points 

We refer ^ for the details of the generation of gauge configurations. We generate 2,500 
trajectories on a 16 3 x 48 lattice for ten combinations of up-down and strange sea quark masses, 
i.e. five m u /s times two m/s. 

We calculate 80 pairs of the lowest-lying eigenmodes on each gauge configuration and store 
them on the disks. These eigenmodes are used to construct the low-mode contribution to the quark 
propagators. The higher-mode contribution is obtained by conventional CG calculation with sig- 
nificantly smaller amount of machine time than the full CG calculation. Those eigenmodes are also 
used to replace the lower-mode contribution in the meson correlation functions by that averaged 
over the source location (low-mode averaging ) 10,01. We extract meson masses from the expo- 
nential decay of the time-separated correlation function of pseudo-scalar operator (P(t)P(O)). The 
decay constant, which is defined by the matrix element of the axial-current operator A^, is obtained 
simultaneously using the PCAC relation d^A^ = 2m q P. 

Throughout the Monte Carlo updates, the global topological charge of the gauge configurations 
is fixed to zero. This is necessary to avoid discontinuous change of the Dirac eigenvalue, which 
is numerically too-expensive. The artifact due to fixing the topology is understood as a finite size 
effect [||] in addition to the conventional one. For the physical size of our lattice L Ri 1.7 fm, the 
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Figure 1: Chiral extrapolation using NNLO full SU(3) ChPT formulae. The left panel contains m^/m u d 
(red circles) and m\ / m sc / (blue squares) while the right contains f K (red circles) and fx (blue squares). Filled 
(open) symbols and solid (dashed) curves indicate the result with the lighter (heavier) strange quark mass. 



finite size effect could be sizable. We calculate both kinds of finite size effect from the analytic 
formulae based on ChPT [Q, g]. In particular, for the effect of the fixed topology, we make use of the 
numerical data of the topological susceptibility determined on the same lattice configurations 

In order to obtain the physical quark mass, we need to renormalize bare quark mass on the 
lattice as m^ ren ' = Z m mq . We obtain the renormalization factor Z m by calculating scalar and 
pseudo-scalar vertex functions in the momentum space in the Landau gauge and applying the 
RI/MOM scheme [10]. In extracting Z m from the vertex functions, we control the contamina- 
tion from the spontaneous chiral symmetry breaking by using the the low-mode contribution to 
the chiral condensate [11]. Using the perturbative matching factor known to 4-loop level and the 
extrapolation to the chiral limit, i.e. m uc i = m s = 0, we obtain the result Z^f s (2 GeV) = 0.806(27). 

In the rest of this article, it is understood that all data points are corrected by the finite size 
effects and quark masses are renormalized. 



3. Fit to NNLO St/ (3) ChPT 

Since we found in the two-flavor calculation that the NNLO ChPT formulae can nicely fit our 
data even in the kaon mass region if one uses the -expansion, we apply the same strategy for our 
2 + 1-flavor analysis. As functions of % K = 2nv\ / '{An f n ) 2 and ^ = 2m\ / '{An f n ) 2 , predictions from 
the SU (3) ChPT are expressed as 

ml/m ud = 2B [l+M^{^ K ;L r 4 ,L r 5 ,L r 6 ,L r 8 )] + a^ ■ + a% ■ Uk + a? ■ & (3.1) 

m\lm sd = 2B [\+M K {^^ K ;L r 4 ,L r 5 ,L r 6 ,L r & )]+af + eg (3.2) 

fn = fo [1 + WA)] + Pi ■ + M ■ Uk + ft ■ & (3.3) 

h = fo [l+F K {^ K ;L r 4l L r 5 )] + j3f • ~ &)+P? • - k), (3.4) 

where m sc j = \{m s + m U( j). When we write the chiral Lagrangean as = + ^4 + H 

with ££ n indicating the contribution of G{p n ). a ( - ' and j8 ( - ' are unknown parameters for the tree- 
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Figure 2: Chiral extrapolation using the prediction from the reduced SU(2) ChPT to NLO. The organization 
is same as Figure 1 . For the fit, the lightest three data points are used for each quantity. 



level contribution from Jzfg- Expressions for the functions M n , M K , F 71 and F K are too involved 
to present here [12]. We note that each function contains NLO contributions from the one-loop 
effect of Jz?2 and tree-level effect of Jz?4 and NNLO ones from the two-loop effect of J2?2 and the 
one-loop effect of These contributions are accompanied by LECs for NLO ChPT, i.e. L[-L r s . 
However, L^L^L^ and L r 7 appear only in the NNLO terms and cannot be determined precisely. We 
introduce values L\ = (0.43 ±0.12) • 10~ 3 , L r 2 = (0.73±0.12) • 10~ 3 , L\ = (-2.53 ±0.37) • 10~ 3 
and L' 7 = (—0.31 ±0.14) • 10~ 3 (defined at ji = 770 MeV) from a phenomenological estimate Jl3| ] 
and determine others L r 4 ,L r 5 ,L r 6 and L\ by a fit. Thus, the chiral extrapolation with (3.1 )— ( 3 .4 ) 
contains 16 fit parameters in total. 

We fit m 2 n /m U( i, m\/m s ^, f % and fg simultaneously taking the correlation within the same sea 
quark mass (m U( t,m s ) into account. By using all data points, reasonable quality of the fit is obtained 
with # 2 /dof = 2.52. In this new study, we determine the lattice scale by the result of f n extrapolated 
to the physical point with the input f n = 130.0 MeV. As a result, we obtain a" 1 = 1.968(39) GeV 
and the pion mass covers the range of 340 MeV < m % < 870 MeV. Figure [j] shows all quantities in 
question as a function of m 2 n . Different symbols represent the pion data (rn\jm uc i and f n ) and the 
kaon data (mjc/m S( i and f£). Filled (open) symbols represent a fixed lighter (heavier) strange quark 
mass, which is accompanied by the solid (dashed) curves. 

Extrapolating the data to the physical point ), which is determined with m K = 

135.0 MeV, m K = 495.0 MeV and f % = 130.0 MeV, we obtain preliminary results m^f (2 GeV) = 
3.64(12) MeV, mp(2 GeV) = 104.5(1.8) MeV, m s /m ud = 28.71(52), f K = 157.3(5.5) MeV and 
Ski 1 fn = 1.210(12), where the errors are statistical only. 



4. Fit to the reduced SU (2) ChPT to NLO 

As a check of the chiral extrapolation with the NNLO ChPT in the previous section, we also 
study a different fit ansatz. It is also possible to carry out the extrapolation to the physical point 
£jf by paying attention only to the dependence on the up-down quark mass, or the pion mass. 
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Integrating out the strange quark as a static heavy quark in ChPT, one obtains an effective theory 
which respects a reduced SU (2) symmetry [14, 1^, [if]]. At NLO, the chiral expansion reads 



m\lm ud = 2fi(l + ±^ln^)+c 3 ^, (4.1) 

m\jm sd = wM+c^Zx, (4.2) 

fn = /(l-^ln^) + c 4 ^, (4.3) 

h = f [K) {l-l^U)+cf ] ^ (4-4) 

where we have LECs B( K \ f( K \ cf ] and cf ] in addition to the LECs for 51/(2) ChPT, i.e. /, fi, 
Z3 and I4. In the present case, all LECs depend on strange quark mass. With the lightest three m u d 
points, which are in the valid region of this framework, i.e. m llt j <C m s for each fixed value of m s , 
we carry out the correlated fit for the quantities sharing the same mass point (m U( j,m s ). Figure || 
shows the fit curves obtained in this way. 

The fit results for each fixed m s are extrapolated to the physical strange quark mass m( phys ' , 
which is determined by solving m K /m s \ ,(phys) = (495.0MeV) 2 /m v . In the light panel of Figure 0, 

Sir 1-1 

we compare physical results for m uc i, m s , fx and fx/ fn from the full NNLO SU (3) ChPT (circles), 
and from the NLO reduced SU (2) ChPT (squares from our analysis and diamonds from the similar 
analysis by RBC and UKQCD Collaborations Jl5|]). The reasonable agreement among different 
fitting prescriptions provides a good consistency check of the analysis. 

It is also interesting to compare SU (2) LECs obtained in this work with our previous work with 
Nf = 2. In the right panel of Figure ||, results of /, £ = Bf 2 /2, I3 and I4 from the Nf = 2+1 sim- 



ulation (circles for the present work and triangles for [|15fl) are compared with the results obtained 
in [§] (squares from the NNLO fit and diamonds from the NLO fit). The reasonable agreements 
observed for each quantity imply that the reduced SU (2) ChPT intermediates full theories with 
Nf = 2 and 2+1. 



5. Summary 

We studied the chiral property of ChPT by comparing the analytic prediction with the lattice 
data obtained in the Nf = 2 + 1 dynamical simulation with the overlap fermions. We fitted the 
data to the ChPT prediction at NNLO for the first time. It is the only way to describe both pion 
and kaon data on an uniform basis and to study the convergence property of the SU (3) ChPT. 
The validity of the extrapolation to the physical mass point is checked with the results from the 
fit with the reduced SU (2) ChPT. However, in order to discuss the convergence property of ChPT 
as in the case of Nf = 2, we need to determine individual LECs of the SU (3) ChPT with good 
accuracy. With the data points obtained for two different strange quark masses, we have a limited 



constraint on the ^ dependence hence large errors for LECs. For example, in (3.3) or (3.4), fo 
cannot be determined unambiguously unless the contribution of the terms with L\ and L r 5 are well 
fixed. Moreover, from the phenomenological side, it is advantageous to determine LECs from this 
work because the results can be used as inputs in the calculation of different quantities including 
Bk and form factors. For these motivations, we are generating more data points with m U( j = m s 
to obtain the SU (3) LECs with high accuracy. 
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Figure 3: Left: comparison of physical values for full QCD obtained by different chiral extrapolations. 
Right: comparison of the SU (2) LECs between different dynamical flavor Nf. All indicated errors are only 
statistical. 



Another work in progress is the calculation on a larger volume for a direct check of finite size 
effect (FSE). We are generating the data for the two lightest ray's and the lighter m s on a 24 3 x 48 
lattice with the same coupling constant as the present work. We are also planning to obtain the 
lattice scale from 12-baryon mass on the larger volume to compare with the current value from the 
/„ = 130 MeV input. 

Numerical simulations are performed on Hitachi SRI 1000 and IBM System Blue Gene Solu- 
tion at High Energy Accelerator Research Organization (KEK) under a support of its Large Scale 
Simulation Program (Nos. 08-05 and 09-05 ). The work of HF was supported by the Global COE 
program of Nagoya University "QFPU" from JSPS and MEXT of Japan. This work is supported 
in part by the Grant-in-Aid of the Ministry of Education (Nos. 19540286, 19740121, 19740160, 
20105001, 20105002, 20105003, 20105005, 20340047, 21105508, 21674002) and the National 
Science Council of Taiwan (Nos. NSC96-2112-M-002-020-MY3, NSC96-2112-M-001-017-MY3, 
NSC98-2119-M-002-001), and NTU-CQSE (Nos. 98R0066-65, 98R0066-69). 
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